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Abstract. We study entropy production and fluctuation relations in the restricted 
solid-on-solid growth model, which is a microscopic realization of the KPZ equation. 
Solving the one dimensional model exactly on a particular line of the phase diagram 
we demonstrate that entropy production quantifies the distance from equilibrium. 
Moreover, as an example of a physically relevant current different from the entropy, we 
study the symmetry of the large deviation function associated with the interface height. 
In a special case of a system of length L = 4 we find that the probability distribution of 
the variation of height has a symmetric large deviation function, displaying a symmetry 
different from the Gallavotti-Cohen symmetry. 



1. Introduction 

The theoretical understanding of out of equilibrium systems is one of the main goals in 
statistical physics. Nonequilibrium phenomena display very rich behavior, with several 
features not observed in equilibrium, as for example phase transitions in one-dimensional 
systems with short-range interactions. Recently great theoretical progress has been 
made concerning systems in contact with two reservoirs [TH5]. Nevertheless, for general 
nonequilibrium systems very little is known. 

For some continuous time Markov jump processes an (average) entropy production 
can be defined [SIE]- This quantity is zero if detailed balance is fulfilled, which 
corresponds to an equilibrium stationary sate, and it is positive in the case of a 
nonequilibrium stationary state. Therefore, a positive entropy production in the 
stationary state is a signature of nonequilibrium. In spite of the generality of this 
statement the relation between entropy production and nonequilibrium stationary states 
is still poorly understood [H]. Recent studies in this direction concern the relation 
between entropy production and nonequilibrium phase transitions, where the entropy 
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production as a function of the control parameter has been observed to peak near the 
critical point p-[TT]. 

In this paper we will consider a microscopic surface growth model in the Kardar- 
Parisi-Zhang [12] (KPZ) universality class. The KPZ equation reads 

dth{x, t) = v + z/V^/i(x, t) + (A/2)(V/i(x, t)f + r/(x, t), (1) 

where x denotes the position on a (i— dimensional interface (we will be dealing only 
with the d = 1 case), h{x,t) represents the height of the interface and T]{x,t) is 
a Gaussian white noise. The Laplacian term is related to surface tension, v is the 
velocity of the interface at zero slope and the non-linear term is the lowest order term 
that breaks the up-down symmetry {h{x,t) — )■ —h{x,t)) [T31II1]. Assuming that a 
nonequilibrium stationary state of a KPZ interface is characterized by v and A, we 
study how the entropy production varies with these parameters. We shall argue that 
this is a reasonable assumption and it allows us to investigate how entropy production 
characterizes a nonequilibrium stationary state. We will work with the restricted solid 
on solid (RSOS) growth model, because this model has a rich phase diagram and we 
can perform analytical calculations (exact and approximative). 

The positivity of the average entropy production is a statement analogous to the 
second law of thermodynamics. More generally, it is possible to associate a fluctuating 
entropy with a single stochastic path [7], and the probability distribution of this 
fluctuating entropy is constrained through the fluctuation relation [T3H21], which is 
a stronger statement than the positivity of the average entropy production. The 
fluctuation relation implies the existence of a symmetry of the large deviation function 
associated with the probability distribution of entropy which is known as the Gallavotti- 
Cohen (GC) symmetry [IZIIIHIES] . 

In general, unlike the entropy production, other time-integrated currents are not 
expected to have symmetric large deviation function [12]. In a growing interface, the 
physically relevant time-integrated current is the variation of the interface height. The 
main focus of this work is to analyze the symmetry of the large deviation function 
associated with the variation of height using the RSOS model. We find that while 
indeed no symmetry of the large deviation function exists in the general case, the large 
deviation function related to the variation of height is symmetric in the case of a four- 
site system. This is in spite of the fact that entropy production and height are not 
proportional to each other in the long time limit. Unlike the GC symmetry which 
results from a certain relation between the rate of each microscopic trajectory and its 
time reversed one, this symmetry is a result of a similar relation between the rates 
associated with certain groups of trajectories. 

The paper is organized as follows. In Sec. 2 we define entropy and briefly discuss 
its relation with nonequilibrium stationary states. In Sec. 3 we define the RSOS model 
and calculate the average entropy production analytically, analyzing how it depends 
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on the interface velocity and A. The fluctuation relation and the GC symmetry are 
introduced in Sec. 4. In Sec. 5 we investigate fluctuations of the interface height. We 
show that if a second time-integrated current is considered a fluctuation relation for the 
joint probability distribution of two currents, with one of the currents being the height, 
can be found and we also show that for the four-site system the large deviation function 
related to height is symmetric, displaying a symmetry different from GC symmetry. We 
conclude in Sec. 6. 



2. Entropy production in nonequilibrium stationary states 



In the following discussion we consider continuous time Markov processes that obey the 
master equation 

jP{C, t)=J2 [wc'^cPiC, t) - wc^c'P{C, t)], (2) 

where P(C, t) is the probability of being in state C at time t and wc^c is the transition 
rate from state C to C . For such stochastic processes three different fluctuating 
entropies can be defined, namely [7] 

(a) the configurational entropy of the system Sgysit) = — In P(C,t), 

(b) the entropy of the medium (=environment) Sm{t) which changes 
instantaneously by \n{wc-^c' /'Wc'-^c) during a transition from C to C. 

(c) and the total entropy Stotit) = Ssysit) + Smit). 

Note that the definition of the environmental entropy requires the rates to satisfy the 
condition that if wc^c > then the reverse rate wc^c has to be positive as well. 

A stochastic path is a sequence of transitions 

Co — Ci — 7- C2 — )■ . . . Cn , 

taking place at random instances of time Tq < ti < t2 < ■ ■ ■ < t^ < T according to 
the specific rates during a fixe time interval AT = T — Tq. Such a path of transitions 
changes the aforementioned entropies by 

A^,,, = In P(Co, to) - In P{Cn, t^) , (3) 

N 

AS^ =^ln^^^^^^, (4) 

and AStot = ASsys + ASm- While these entropies vary discontinuously for a particular 
stochastic path, their expectation values vary smoothly according to 

Ssysit) := J^iSsysit)) = Y.^^ {C , t)wc^c' , (5) 

C,C' ^ ' ^ 
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Figure 1. Deposition and evaporation rates for the RSOS model. 

Stot{t) 

where (. . .) denotes the ensemble average. Since the last equation can be rewritten as 
stotit) = i5^(ln[P(C,t)ii;c^c']-ln[P(C",t)ti;c'^c]) (^P{C,t)wc^c'-P{C' ,t)wc'^c){8) 

C,C" 

and both brackets have always the same sign, the average total entropy can only increase 
and will stay constant if detailed balance is satisfied, in accordance with the second law 
of thermodynamics. Moreover, if the system is in a non-equilibrium stationary state, the 
average internal entropy (Ssysit)) = — P(C) In P(C) is constant so that Stot = Sm- 
Since the stationary master equation 

= wc'~,cP{C') - J2 wc'-.cPiC) (9) 
c c 

is invariant under the replacement wc-^c' ^ "W^c-i-c" + ^{C, C')/ P{C) with an arbitrary 
symmetric function A{C, C) = A{C', C) while Eq. (|8]) is not invariant under this 
operation, the same stationary state may be generated by different dynamical rules 
with a different entropy production rate (see [8]). 



dt 

d_ 

dt 



(Srait)) = y2ln^^^^P{C,t)Wc-.C 



C,C' 



(6) 



P\C', t)wc'-,c 



3. Entropy production in a solid-on-solid growth model 
3.1. Definition of the model 

As a microscopic realization of the KPZ equation let us consider a restricted solid- 
on-solid (RSOS) growth model on a one-dimensional lattice with L sites and periodic 
boundary conditions. The configuration C of the interface is characterized by height 
variables hi eTj attached to the lattice sites i which obey the restriction 

\h,-h,±i\<l. (10) 

The model evolves random-sequentially by deposition and evaporation of particles with 
rates qa,b,c and Pa,b,c according to the rules shown in Fig. [1] Because of (fTOjl . it is 
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convenient to describe the configurations of tlie interface in terms of charges 

ai ■= hi+i-hi e {0,±1} (11) 
which evolve random- sequentially according to the rules 



qa qt qc qc 

; nn ^ -i— n+ , ' +n -0:;zz 



(12) 



Pa Pb Pc Pc 



Since deposition and evaporation correspond to negative and positive charge 
displacements, the average interface velocity v = ji{hi) is given by 



v = {qb-Pa)m +go(-+) + gc((0+) + (-0)) -pe((+0) + (0-)), (13) 

where (ajCTj+i) is the probability to find two specific charges at neighboring sites. 
Likewise the entropy exported to the environment is given by 



U In ^ - p Jn (00) + gjn^(-+) - p,ln*(+-) 

V Ph Pa^ Pa Pb 

+ q,\n^ ((0+) + (-0)) - pM- ((+0) + (0-)) . (14) 



The interface slope m is defined by 

L 

L 



m = j^ai. (15) 

Since we consider periodic boundary conditions we have m = and, therefore, the 
number of positive charges is equal to the number of negative charges. One way 
of introducing a non-zero interface slope in the RSOS model is to change the RSOS 
restriction between sites L and 1 so that hi — hi = H, H ±1. In this case the interface 
slope is given by m = H/L. It can be shown that the A term in the KPZ equation can 
be defined for a microscopic model. It is related to how the velocity depends on the 
interface slope in the following way |13] . 



v{m) = V + {X/2)m^. (16) 

We assume that a nonequilibrium stationary state of the RSOS model is 
characterized by v and A, this means that it is characterized by the interface velocity 
and by how the velocity depends on the interface slope. In the following we calculate the 
entropy production, the velocity and the parameter A defined above. We are interested 
in investigating how the entropy production is a function of these two quantities. 
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Figure 2. Left: Phase diagram of the RSOS model obtamed from numerical 
simulations and the line ([T5| where the model is exactly solvable. Right: Mean field 
phase diagram (see text). 



3.2. Entropy production: Exact results 

As shown in Appendix A, the RSOS model can be solved exactly in the thermodynamic 
limit L — 7- oo in a certain parameter subspace |25]. To simplify the analysis, we will 
from now on restrict ourselves to a two-dimensional section of the parameter space 



qa = qb = qc = q Pa=p Pb = Pc = ^ (17) 

controlled by two parameters q and p. In this section the exactly solvable subspace 
corresponds to the line 

As shown in the phase diagram in the left panel of Fig. [21 this line intersects with the 
lines for = and A = at the point p = g = 1, where the model can be shown to 
obey detailed balance [26]. Note that the lines w = and A = intersecting at the 
detailed balance point is in agreement with the assumption that v and A characterize 
the nonequilibrium stationary state of the RSOS model. 

Along the exactly solvable line, the density of charges (cf. Eq. (]A.12p ) is given by 

2 . ^ 

P = , 19 

2+ v/(2g-l)/g 

As shown in appendix B, the two-point function in the charge basis factorizes in the 
exactly solvable line, so that the velocity w, the nonlinearity coefficient A, and the entropy 
production rate Sm are given by 



(g-l)p, (20) 
4g(l - p) 



A =-(g-l)-— ^ ^ i + 6g-2p(2g+l) , (21) 
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Figure 3. v{q),X{q) and Sm{q) along the exactly solvable line. 

Sm. = iq-l) (y In + |(2 - p) In g) . (22) 

In Fig [3] we plot the three quantities as a function of q. We see that the entropy 
production is always positive and it grows as the distance from the equilibrium point 
q = I increases. The velocity and A have opposite signs and as their absolute values 
increase the entropy production increases. This is in agreement with the idea that the 
entropy production gives a measure of how far the system is from equilibrium. 



3. 3. Entropy production: Approximate results 

The two-point function in the charge basis factorizes only on the line p = q{2—q) / (2g— 1). 
In order to access the whole phase diagram we now assume that it factorizes for any 
p, q, which is analogous to a pair-mean field approach in the height representation 
j. Details of the calculation are presented in Appendix B. Within this mean-field 



approximation, from fIB.Sj) . the density of charges in the stationary state (at zero slope) 
is given by 

p = 2/(2 + ((g + p)/(g + l))-^/^). (23) 
Repeating the above calculations we find that the interface velocity vanishes for 

P= -q + q'^ + q^- (24) 
and that the coefficient of the non-linear term in the KPZ equation vanishes for 



p = (-3 - g + + ± V5 + lOg - 9g2 - 20g3 + 3g4 + lOg^ + g6)/8 , (25) 



where the sign in front of the square root is positive for p > 1 and negative ioi p < 1. 
As shown in the right panel of Fig. [21 these lines approximate the true lines, cross at 
the equilibrium point, and have the same type of curvature. In addition, the figure 
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shows the hne where the entropy production obtained form equations f l23|) and f lB.Sp . 
is minimal. This hne does not coincide with the hne w = although it is very close 
to it. Similar situations have been observed in other models [OHTTj. where the entropy 
production was found to peak near the critical point. In the present case, this happens 
because although the entropy production strongly depends on w, it also depends on A. 

In principle Sm can be considered functions of v and A. However, we could find 
analytical expressions only for w = 0, where p = —q + q"^ + q^. On this line, from (125]) 
and (1B.7|) we find 

X = {q-l)/q. (26) 

For large q we have A — ?■ 1 and at g = (v^ — l)/2), which corresponds to p = 0, 
A = -(VS - l)/2 ^ -0.618. From (ES]) we have 

--^■"(^)- P^) 

where the above relation is valid for — (-\/5 — l)/2 < A < 1. As expected, the entropy 
production is zero at A = and grows as the absolute value of A grows. Note that it 
diverges for q — )■ (v^ — l)/2) and g — j- oo. Again the results are in agreement with 
the idea that entropy production gives a quantitative measure of the distance from 
equilibrium. 



4. The Gallavotti-Cohen symmetry 

The fluctuation relation states that 

where P{ASm) denotes probability of the variation of the medium entropy during time 
AT and the above relation is valid for AT — oo. As we show below it leads to the 
GC symmetry, which is a symmetry in the large deviation deviation function associated 
with P{ASm)- In order to do that in a more general framework we first introduce 
time-integrated currents and the large deviation function associated with them. 

Time-integrated currents are a functional of the stochastic path, which changes 
only when the underlying process jumps to a different configuration. They are defined 
by an increment Oq^c' that depends on the transition rates and is antisymmetric, i.e. 
Oc^c = —(^c'^c- As before we consider the path Co — Ci — t- ... — t- Cn starting at time 
To and finishing at time T, where is the total number of jumps and AT = T — To. A 
time-integrated current takes the generic form 

N 

J = Y,dc,^.^c,. (29) 

fc=i 
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In the case Qc^c = In "'g^c ^ J is just the entropy change ASm Q- Another current 
we are interested in is the variation of height, where 9c^c' is +1 for a deposition and 
— 1 for an evaporation. 

The average of the time integrated current (l29l) is 

pAT 

(j)= / dtY^o c-^c'P{C,t)wc^c' (30) 

° C,C' 

In the long time limit, we have by ergodicity 
J 



AT 

c,c 



J2^c^c'P{C)wc^c', (31) 



where P{C) is the stationary state probability distribution. The large deviation function 
associated with this time-integrated current is related to deviations of the current from 
its average value. The large deviation function is defined by 



P(j)^exp[-ATJ(j)], (32) 



where P{j) denotes the probability that j = and the above relation is valid for 
AT — )■ oo. The GC symmetry is a symmetry of the large deviation function that 
reads [T7| 

/(j) - li-j) = -Ej (33) 

The fluctuation relation for the external medium entropy (128|) corresponds to E = 1 
and was first proved, for jump processes, by Lebowitz and Spohn [19]. The central idea 
is to introduce the deformed time evolution operator 

Note that for s = we get the usual time evolution operator for the master equation. 
Consider the joint probability Pj{C,T) of being at a state C at time T with time- 
integrated current J. The Laplace transform of it is 

P,(C,T) = J2PAC,T)exp{-sJ). (35) 
J 

From the master equation ([2]) we obtain [19] 

4f^PsiC,T) = -Y,Hcc'{s)Ps{C\T). (36) 

Since Hi^s) is a Perron- Frobenius operator, it is possible to relate the generating function 
of the time-integrated current J with its minimum eigenvalue in the following way, 

(exp(-sJ)) ~ exp(-AT/(s)), (37) 
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where J(s) is the minimum eigenvalue of the matrix H[s). The Gartner- Elhs theorem 
states that the large deviation function (1521) is the Legendre transform of /(s), 



/(j) = mm [I{s) - sj) . (38) 

Therefore the GC symmetry (l33l) is equivalent to 

i{s) = i{E-s). (39) 

In the case where the time integrated current is the external medium entropy (jl]), we 
have 



Hcc'{s) = { ^"^^'r.^"^^'-^" , (40) 



Y,c" ^c^c" if C = C 

which clearly has the symmetry H^{s) = H{1 — s). From this condition, follows fl39l) 
with E = 1. 

Let w'^^fj, be some equilibrium transition rates related to an equilibrium 
distribution P'^'^(C), where detailed balance is satisfied. A sufficient condition for a time- 
integrated current corresponding to transition rates wc^c to satisfy the GC symmetry 
is that a quantity E exists, such that these rates take the form [23] 

E 

wc^c' = w^c^c' exp(— 6'c^c')- (41) 

Equivalently to (SI]), we have 

wc^c = P;,\C)wc'^cPe,{C') expiEOc^c') (42) 

Defining Peg as a diagonal matrix with elements Peq{C), it then follows, from the last 
relation, 

H^is) = P-g'HiE - s)Peg, (43) 

which directly gives (l39l) . We remark that time-integrated currents, such that an E 
satisfying ( HTj) can be found, are proportional to the entropy (up to boundary terms) in 
the long time limit (see [l9l[3T]). More precisely, 

k=l ^ ^ ^' i=l 

Therefore, condition ( 14T]) is very restrictive in the sense that very few currents satisfy 
it. 

Nevertheless, if one considers the joint probability of more than one time-integrated 
current a more general relation can be found [19]. Let us take k currents. Similar to 
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fj^Tj) . a sufficient condition for a fluctuation relation including k currents is that the rates 
take the form 

(45) 

We define the large deviation function associated with the joint probability distribution 
of the currents Jj by 

m,J2,...,Jfc) ~exp(-AT/(ji,j2,...,Jfc)). (46) 
Analogously to flMl) we can define the deformed operator 

S- , ._\-wc'^c(iM-Y!l=iStOh^c)'^^C^C' 

tlcc'[.Si,S2,...Sk) - \ ^ -f n n, > (47) 

1^ l^C" ^c^c" It o = o 

and relation fH5|) gives 

H^{si, S2, ...Sk) = P-'H{E, - si, E2 - S2, E, - Sk)Per (48) 
From the last relation we have a multidimensional version of (139|) . which reads 

J(Si, S2, ...Sfc) = T{Ei - Si, ^2 - S2, ^fc - Sfc). (49) 

In terms of the large deviation function we have the following generalized GC symmetry 

k 

^(ji,j2,---,Jfc) - /(-ji, -j2, -jfc) = - ^Eiji. (50) 
5. Fluctuations of the interface height 

Let us consider the RSOS model again. The physically meaningful time-integrated 
current is the variation of height. Along the line p = 1 the height is just the entropy 
multiplying In q. Therefore, for p = 1 the probability distribution of the variation of 
height displays the GC symmetry fl33|) with E = \nq. For p ^ 1 the relation between 
height and entropy becomes more complicated, and we would like to investigate the 
symmetry of the large deviation function associated with the probability distribution 
of height for this case. In other words, we would like to know if a symmetric large 
deviation function (related to height) is a general property of the RSOS model or if it 
holds only for specific microscopic rules {p = 1). In the rest of this section we show 
that if a second time-integrated current is considered, fluctuation relations for the joint 
probability of two time-integrated currents that holds in the whole phase diagram can 
be found. We then proceed to analyze the symmetry of the Legendre transform of large 
deviation function associated with height for small systems [L = 3,4). 
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5.1. Joint probability of two time-integrated currents 

We introduce the time-integrated currents AHq and AHgp which are defined in the 
following way. The increment 6c_^c/ of AHg is —1 for an evaporation taking place with 
rate 1, +1 for the deposition related to the reverse transition, and otherwise. In a 
same way, the increment O'^^^j, of AHgp is —1 for an evaporation taking place with rate 
p, +1 for the deposition related to the reverse transition and otherwise. Note that, in 
the present case f l45p is verified on the form 

wc-^c' = wc'^c-exp (^(Inq) 9^^^, + (^n ^ O^^^^^ . (51) 
Then, relation fl50|l becomes 

HK^ Kp) - H-hq^ -hqp) = - ln(g)/ig " \n{q / p)hqp (52) 

In order to include the height as a current we note that in terms of AHq and AHgp, the 
height variation (multiplied by the system size L), AH, is given by 

AH = AHq + AHqp (53) 

Hence, fl52]) implies in 

I{h, hqp) - I{-h, -hqp) = - ln{q)h + \n{p)hqp 

lih, hq) - I{-h, ~hq) = - \n{q/p)h + \n{p)hq. (54) 

The height appears explicitly in these relations and they are valid in the entire phase 
diagram. Therefore, it is possible to find fiuctuation relations for the height if a second 
time-integrated current is considered. Next, we consider the probability distribution of 
the variation of height, without a second time-integrated current. 



5.2. A symmetry different from the GC symmetry 

We consider the RSOS model for very small system sizes starting with L = 3, where 
the model has 7 states in the charge representation. In order to visualize the network of 
transitions in this space, let us denote by [ai, . . . , a^] all configurations which differ only 
by translation, e.g. [^ — 0] = {^ — 0, OH — , —0+}. As shown in the left panel of Fig. HI 
the transitions for L = 3 form a simple cycle. Although each transition contributes 
differently to AH and ASm. it is clear from the value of the rates (i.e. Fig. H]) that for 
each completed cycle the resulting changes will have the fixed ratio ^^ff- = ln(g) — | ln(p). 
Therefore, in the long-time limit AH will be distributed in the same way as ASm and 
thus trivially displays the GC symmetry. More precisely, it is possible to find a similarity 
transformation of the form (143 p with E = ln(g) — |ln(p). 

For L = 4 the situation is different. As shown in the right panel of Fig. HJ the 
transition network is a loop consisting of two smaller cycles. If the system passes the 
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L=3: [jo^] 

3q 3p 

1 q 

^ 2q ^ 

[+-0] ^ " [+0-] 
2 



L=4: ^[0000]^ 



'[--++] 



^[+-00] < 



,[+0-0] 



'[+-+-] 



[+00-]- 



Figure 4. Network of transitions in the charge representation of the RSOS model 
with L = 3 (left) and L = 4 sites (right). 



whole loop it produces both entropy and height, but if it circulates within the small 
cycles it generates only entropy. This means that entropy and height are no longer 
coupled and are expected to be distributed differently. 

To find out whether the distribution of Aif is still symmetric for L = 4, let us 
follow the procedure described in the previous section. In the charge basis the deformed 
evolution operator H{s) turns out to be a 19 x 19 matrix of which the lowest eigenvalue 
has to be found. The problem can be simplified even further by defining the 6x6 matrix 

T.A^) := ^^ec'(s), (55) 

where k and k' denote the six sets of charge configurations which differ only by 

translations and C is one of the states represented by k'. In the basis {[0000], [ 

++], [+ - 00], [+ - +-], [+0 - 0], [+00-]}, this matrix reads 

\ 

— 1 — g pe'^ qe~^ 

Aqe-' qe"" -l-j9-3g 2e^ 2e' 

Tnn'{,s) = ^ ^ . . . (56) 



4p — Aq 





e' 








qe 





-1-q 


pe'^ 








qe~ 


4ge~'^ 


qe"'^ 


— 1 — p — 3q 


2e' 


2e' 











qe""" 


-2 - 2q 





pe^ 








2qe-' 





-2-2q 


2e' 


4pgS 







2qe-' 


2qe~' 


-2-p 



V 

Since the eigenvector corresponding to the lowest eigenvalue of H{s) is translationally 
invariant, there will be a corresponding eigenvector of T{s) with the same eigenvalue. 
Hence I{s) can be computed by determining the lowest eigenvalue of the matrix given 
above. 

Computing its determinant we find that I{s) vanishes for s = and E' := s = 
^(ln3g^ — ln(2p + p^)) and that the characteristic polynomial and hence the spectrum 
is invariant under the replacement s ^ E' — s. Therefore, the RSOS model with L = 4 
sites has a symmetric large deviation function (associated with height). Moreover, the 
function /(s) for the height current - even when appropriately rescaled - differs from 
the corresponding function for the entropy (see Fig. [5]), showing that entropy and height 
are not proportional to each other for AT — > oo. In Appendix C we present a different 
proof for this symmetry, without using the matrix H{s). 

The GC symmetry has a clear physical interpretation. It is a result of the fact that 
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Figure 5. Numerically determined Legendre transform /(s) of the large deviation 
function of the entropy (dashed line) rescaled by s -M^p+p ) ^ g^j^j ^^^^ ^j-^g 

height (solid line) in the RSOS model with p = 0.1 and q = 2. 

the exponential of the variation of the entropy for a given stochastic path is the weight 
of the path divided by the weight of the time- reversed path (see [HIES]). Clearly, this 
is not the case of the height in the RSOS model, therefore the symmetry we found for 
the L = 4 case is different from the GC symmetry. In Appendix C we show that the 
symmetry is a result of some subtle arrangement of the stochastic trajectories that leads 
to a variation of height AH. Hence, we have shown that currents different from the 
entropy may also have a symmetric large deviation function. 

Unfortunately, this observation does not extend to L > 4, where a small asymmetry 
in I{s) can be observed (we calculated I{s) numerically, directly form the matrix H{s), 
for L = 5, 6, 7). This may be caused by the increasingly complexity of the network of 
transitions (see Appendix C). 

6. Conclusions 

In the present paper we studied the entropy production and fluctuation relations for 
the particular solid-on-solid growth model which belongs to the KPZ universality class. 
On a particular line of the phase diagram, where the model can be solved exactly, we 
calculated the entropy production, the velocity, and KPZ-parameter A. The results are 
compatible with the idea that |A| and can be used as a measure for the distance from 
equilibrium. 

Extending these exact results to the entire phase diagram within a mean field 
approximation, we obtained analytical expressions for the lines v = and A = which 
are close to the true lines and exhibit the same curvature. Furthermore, we calculated 
the entropy production as a function of A along the line v = 0. It was found that Sm 
grows with |A|, again in agreement with the idea that the entropy production gives a 
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measure of the distance form criticality. 

As we saw, along the hne p = 1 height is proportional to entropy and therefore 
its probability distribution displays the GC symmetry. We showed that if a second 
time-integrated current is considered a symmetry in the joint probability distribution of 
two time integrated currents valid in the whole phase diagram can be found. Moreover, 
we explored the symmetry of the Legendre transform of the large deviation function 
associated with the variation of height for the L = 3, 4 cases. For L = 3 the height 
turned out to be asymptotically proportional to the entropy, hence it trivially displays 
a GC symmetry. A different result emerged for L = 4: entropy and height are not 
proportional in the long time limit, nevertheless the large deviation function related to 
height is still symmetric. As we showed this a new symmetry different from the GC 
symmetry. 

We would like to point out that the fluctuation relation for the entropy production 
is very general, holding for any Markov process with reversible transitions rates. 
Contrarily, it is valid for very specific currents, that is, currents that become proportional 
to the entropy in the long time limit. Here we found an example of a current different 
from the entropy with a symmetric large deviation function. It would be interesting 
to investigate what currents in Markov processes, different from the entropy, have a 
symmetric large deviation function and what may be the physical origin of the symmetry. 
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Appendix A. Exact calculation of the stationary state of the RSOS model 

Here we provide more details on how the stationary state of the RSOS model expressed 
in the charge basis can be computed exactly in a certain subspace of the parameters by 
using the method introduced by Neergaard and den Nijs [2S]. Let us assume that the 
probability to find the stationary system is proportional to 



fM,E = x''y^ , 



(A.l) 



where x < 1 and y < 1 are certain numbers depending on the rates and 





(A.2) 



3 
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are the total number of the charges (divided by two) and their electrostatic energy, 
respectively. Inserting this ansatz into the master equation yields 

= n+_[gb/M-i,£;+i - PhfuA + n-+[pafM~i,E~i - qafhiA 

+ nQQlqafhl+l^E+l + PbfM+l,E-l — {<lb + Pa)fM,E\ 

+ (no- + n+o) [qcfM,E-i - PcfuA (A. 3) 

+ (ra_o + nQ+)[pcfM,E+l - QcfkLE] 

where n^^ji is the probability of finding a pair aa' of nearest neighbors sites in a 
configuration with M charges and electrostatic energy E. Using the relation 



I M/L ifcr = ±l 

= 1 1-2M/L ifa = ^^-^^ 



one obtains the constraint 

no- + n+o + 2n+_ = n_o + no+ + 2n_+ . (A. 5) 

meaning that there are four independent variables in the master equation (lA.4p . namely 
noo, n+_, n_+ and (no_ + n+o + 2n+_), and that the corresponding coefficients have to 
vanish. The coefficient of (no_ + n+o + 2n+_) vanishes if 

{qc-Pcy^^){l-y) = ^. (A.6) 

This equation has two solutions. The first one, y = Qc/pc, combined with the other 
constraints, reduces to the condition of detailed balance (see The second solution 

y = 1, on which we will focus in the following, is much more interesting because it 
extends to the nonequilibrium regime of the model. In this case the configurational 
probability distribution becomes independent of the electrostatic energy E. 

For this solution the coefficient of noo gives 

x=*^. (A.7) 

Qa+Pb 

The two remaining coefficients of n_| and n |_ lead to the same relation, namely 

qc-Pc = —, ^ ^. A.8 

Therefore, in the region of the phase space where the constraint flA.Sp is satisfied the 
probability of a configuration with M charges in system of size L in the stationary state 
is exactly given by 



where 



Pm = Zll.x^ (A.9) 



= 5Z M\M\{L-2M)f'' ^^-^^^ 

M=0 ^ ' 
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is the partition function and x is given by flA.7|) . In finite systems neighboring sites are 



still correlated because of charge conservation. For example, in a system with 2 sites the 
only possible configurations are 00, +-, and -+ and thus their probability distribution 
does not factorize. However, in the limit L — >■ oo this conservation law will be less and 
less important, meaning that the state factorizes asymptotically: 

{aa') = {a){a')+oil/L), (A.ll) 

Moreover, in the large L limit the sum in (lA.lOl) is dominated by its maximum and the 
density of charges p = 2M/L is given by 

(A.12) 



2 + X-1/2- 

Appendix B. Mean field calculation of A and Sm 

To compute the value of the coefficient A in the KPZ equation one has to consider a 
tilted interface in the stationary state and to investigate how the interface velocity v 
varies with the slope m. According to equation ( 1T6|) . A is then given by 



A = --^v(m) (B.l) 
m dm 

where v{m) is given by Eq. f irS]) . 

To compute v{m) we assume that CTjCTj+i factorizes. Of course this is only correct 
for m = in the parameter subspace constrained by ( ]A.8|) . However, it is reasonable 
to expect that the stationary state still factorizes if m is sufficiently small. With this 
assumption we have 

, , p(m) + m , , n(m) — m , , , , x 

(+) = ^ — ' (-) = ^ — ' (o) = i-pM (B.2) 

with some unkown function p{m) for the charge density, so that the interface 
velocity (IT^ is given by 

v{m) = {qb - Pa){l - p{m)f + ^(g^ - Pb){p{mf - m?) 

+ {<lc~Pc){l- p{m))p{m); (B.3) 
To compute A we need to know the derivative of p{m). Since 

j^p{t) = {qb + Pa)m-qa{-+)-Pb{+-) (B.4) 
vanishes in the stationary state, we have 

{qb + Pa){l- p{m)f -h^qa+Pb){p{mf -m") = 0. (B.5) 
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Differentiating with respect to m we obtain the expression 



m —p{m) = p + 4 ■ . 

dm qa + Pb 

which tends to in the hmit m — )■ 0. Inserting this result into (IB.ip yields 



(B.6) 



A 



X 



-1/2 



4 



■[{qa - Pb)p -4:{qb-Pa){l- p) + 2(gc - Pc)(l - 2p)] 



{qa - Pb) 



(B.7) 



One can obtain the same expression for A by differentiating the deterministic part of 
the KPZ equation for the RSOS model [25]. Similarly, we can compute the entropy 
production rate 



{qa In — - In — )^ + {qt In — - p„ In — )(1 - pY + [q^ - Pc) In —p{l - p).(B. 

Pa Pb 4 Pb Pa Pc 



Thus we have obtained the quantities v, Sm and A as functions of the density of charges 
p in the stationary state, assuming factorization of {era'). These results are exact in 
the parameter subspace constrained by f lA.Sp in the limit L — )■ oo, where the density of 
charges is given by flA.12p . 

Appendix C. The special symmetry for L = 4 



In this Appendix we show that the large deviation function associated with the height 
variation I{h) = — hmAr-i-oo lnP(/i), where h = AH/ AT, has the symmetry 



1(h) - li-h) = -h- In ^ , 



(C.l) 



for the RSOS model with L = 4. 



We start by considering the embedding Markov chain (discrete time) of the 
Markovian process. The transition matrix reads 



A 



c 



(C.2) 



where Ac = 'Ylcy^c'^c^c' ■ Within the basis used in (156|) . the transition matrix takes 
the explicit form 



/ 



_2_ 

<?+? 





p 






1 










3g+p+l 





P 








3g+p+l 


q 





1 


1 



















3g+p+l 





2q 








■iq+p+1 


1 





q 


q 


g+1 


9+1 


9+1 



2q+p+2 

q 

2q+p+2 



p 

2q+p+2 
2 

2q+p+2 





(C.3) 
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In this discrete time case we are interested in the probabihty of having a variation of 
height AiJ with a fixed number of jumps A^, which we denote by P( Ai7| A^). Our aim 
is to obtain all possible paths that lead to a variation of height AH within jumps. 
Let us assume that Aif and N are large, such that only periodic sequence of jumps 
contribute, while any other sequence of jumps will just generate a boundary term that 
is irrelevant for large Aif and N . We also consider an initial state [H — 00], which does 
not impose any restriction since we assume Aif and large. 

Among the possible periodic sequence of jumps, there are four cycles that increase 
the height. They are: 

(i) [+ - 00] ^ [+0 - 0] ^ [+00-] ^ [0000] ^ [+ - 00] 

(ii) [+ - 00] ^ [+ - +-] [+00-] ^ [0000] ^ [+ - 00] 

(iii) [+ - 00] ^ [+0 - 0] ^ [+00-] ^ [- - ++] ^ [+ - 00] 

(iv) [+ - 00] ^ [+ - +-] ^ [+00-] ^ [- - ++] ^ [+ - 00] 

These cycles take 4 jumps each and increase height by 4. Defining A = (3g + p + 
l){q + l){q + p){2 + 2q+p) and B = {3q + p + l){q + 1)^(2 + 2q + p), the probability 
of cycle 1 is Ml = 2q^/A, of cycle 2 is M2 = g^/A, of cycle 3 is M3 = 2q^/B and 
of cycle 4 is M4 = q^/B. The reversed cycles decrease height by 4 and they have 
probabilities Mi = p^ /A for cycles 1, M2 = 2p/A for cycles 2, M3 = p'^/B for cycles 
3 and M4 = 2p/B for cycle 4. We denote the number of these cycles in some path 
C = Co — Ci — 7- ... — )■ Cat by mj and the number of reversed cycles by rhi, with 
i = 1,2,3,4. 

Besides these four cycles (and their four reversed ones) there are many other 
periodic sequence of jumps, however they do not change height. For example, it 
is possible to jump from [+0 — 0] to [+00—] and then jump back. This sequence 
takes two jumps and it does not change height. Another example, is the cycle 

[+ - 00] -> [+0 - 0] ^ [+00-] [H \ — ] [+ - 00], which takes four jumps and 

also does not change height. We denote by A the set of independent periodic sequences 
of jumps that do not change height. Each element of the set, i G A, is characterized 
by its probability Pi and the number of jumps it takes 77,. The random number Ui gives 
the number of time the sequence z G A appear in the path {C}. 

The probability of the variation of height AH with a fixed number of jumps A^ can 
be written in the form 

P{AH\N) = J:[^^W[{C}], (C.4) 

where W^[{C}] is the probability of the path {C} and the sum Yl'{c} paths 
with A^ jumps and variation of height AH. By commutativity of the measure of the 
stochastic path, one can order the terms contributing to the probability of a path such 
that a given trajectory looks like a sequence of elements of A, followed by a sequence 
of cycle 1, 2, 3 or 4 (or their reversed). If we have mi + rn2 + fn^ + rn^ = j, then 
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mi + 7712 + ms + m4 = AiJ/4 + j and X]iGA'^*'7«~ ^ ~ ~ ^3- Therefore, the sum 
flC.4p becomes 

iN-AH)/8 I / 

i=o Vie A 

(Ml + M2 + Ms + M4)'^^/^+^' (Mi + M2 + Ma + M^y . (C.5) 

where the term (Mi + M2 + M3 + M4)^-'^/^"''-' comes form the sum over all possible 
mi,m2,m3,m4 subjected to the constraint mi +7712+^3+7714 = AH/A+j and, similarly, 
the term (Mi + M2 + M3 + M4)-' comes form the sum over all possible 77ii, frvz, ms, 7774 
subjected to the constraint ml + m2 + fn^ + 772^ = j. Analogously, for P{—AH\ N) we 
obtain 

(N-AH)/8 / / 

J=o \{vi} VieA 

(Ml + M2 + Ms + M4)^ i^Mi + M2 + Ms + M4j . (C.6) 

Therefore, the probability of the variation of height with a fixed number of jumps has 
the symmetry 

P( AH\ N) _ fMi + M2 + Ms + M4^^ "^^Z' ,^ 




P{ -AH\ N) VMi + M2 + Ms + M4 

Note that right hand side of the last equation does not depend on the number of jumps 
A^. By applying Bayes formula, we have 

P{h) _ j:N=iTh^Pi^H\N)P{N) _ fMi + M2 + M, + M,\^'' 



P{-h) EN=iTh^ P{-m N)P{N) VMi + M2 + Ms + M 

where h = AH/ AT, P{h) is the probability in the continuous time process with fixed 
AT and P{N) is the probability of having jumps. 

Finally, since 

M1 + M2 + M3 + M4 = Q + -^^ and M1 + M2 + M3 + M4 = (p' + 2p)(^ + -^), (C.9) 

from (1C.8|) we obtain the symmetry f lC.ip . 

An important point in the calculation we performed here is that the chosen initial 
state [H — 00] has the property that any cycle that increases height has to go through 
it (we could also choose [+00—]). Without this property it would not be possible to 
write down equations (IC.Sp and (lC.6p . from which we derive the symmetry. As we 
mentioned, the symmetry is lost for L = 5 where the large deviation function displays 
a small asymmetry. By drawing the network of states for the L = 5 case, which has 10 
states, we observed that there is no state such that a cycle that increases height has to 
go trough it. 
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